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ABSTRACT 


A  model  for  the  prediction  of  the  elastic  response  of  composite  materials 
over  wide  ranges  of  concentration  is  presented.  The  method  is  based  on  the 
mathematical  analogy  between  the  motion  of  particles  suspended  in  viscous  media 
and  the  elastic  deformation  of  solids.  The  system  is  locally  described  by 
linearized  forms  of  the  elastic  moduli.  The  validity  of  these  relations  is  then 
extended  to  all  concentrati ons  under  the  assumption  that  any  new  portion  of 
filler  "sees"  the  existing  composite  as  a  non-interacting  homogeneous  matrix. 

The  method  predicts  the  behavior  of  composites  with  solid  spherical  inclusions 
and  foams  over  wide  ranges  of  concentration.  The  model,  free  from  adjustable 
parameters,  shows  excellent  agreement  with  existing  experimental  data.  The 
extension  of  the  method  to  other  inclusion  geometries  is  straightforward. 


INTRODUCTION 


In  the  search  of  a  mathematical  model  for  the  mechanical  behavior  of  hetero¬ 
geneous  systems,  a  complete  analogy  is  found  between  the  basic  equations  for 
elasticity  theory  and  those  for  the  motion  of  suspended  particles  in  viscous 
media.  The  analogy  can  be  easily  proved  by  writing  the  equations  of  motion  for 
an  incompressible  material  in  terms  of  displacements: 


-p  .  +  Gu  •  .  .  +  F . 
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where  G  is  the  shear  modulus,  u  •  is  the  displacement,  F-j  is  the  body  force,  p 
is  the  reactive  pressure,  and  p  is  the  density. 

For  comparison,  the  Navi er-Stokes  equations  for  an  incompressible  Newtonian 

fluid  are: 


’  0  i  n  +  vjh 


where  n  is  the  viscosity,  v  the  velocity,  and  F  the  body  force  per  unit  mass. 

Equations  1  and  2  have  a  term  by  term  equivalence,  with  the  exception  of  the 

nonlinear  component  v-v.  .  representing  a  convective  acceleration.  Under 

J  '  » J 

creeping  flow  conditions,  acceleration  terms  become  irrelevant  and  a  complete 
analogy  can  then  be  established.  Given  a  certain  geometry,  the  theoretical  equations 
for  a  viscous  flow  problem  and  those  for  the  elastic  deformation  of  a  solid  are 
of  the  same  form: 

shear  rate  :  viscosity 
shear  strain  :  shear  modulus 


For  a  suspension  of  rigid  particles  in  an  incompressible  matrix,  there  is  a 
simple  relationship  between  relative  viscosities  and  relative  shear  moduli 


-2- 


(relative  to  the  unfilled  matrix): 


n  G 
m  m 


Thus,  if  a  theory  for  the  viscosity  of  a  filled  system  or  suspension  is 
available,  it  can  be  used  to  estimate  the  shear  modulus  of  an  elastic  system. 

A  second  analogy  can  be  also  found  in  relation  to  the  methods  used  in  the 
solution  of  both  problems,  as  indicated  in  Figures  1  and  2.  In  very  general 
terms,  it  is  found  that  there  is  an  important  group  of  contri buti ons  based  on 
the  establ i shment  of  an  idealized  geometry  and  packing  arrangement  (top  of 
Figures  1  and  2).  A  second  group  is  based  on  the  determination  of  bounds,  which 
avoids  (or  reduces)  the  idealization  problem  by  use  of  energy  and  variational 
techniques  (bottom  of  Figures  1  and  2).  Typical  limitations  of  geometric  models 
and  oound  methods  gave  origin  to  a  third  group  of  contributions  containing  ele¬ 
ments  of  both  techniques  (middle  of  Figures  1  and  2). 

With  regard  to  the  problem  of  predicting  the  elastic  performance  of  com¬ 
posite  materials,  (Figure  2),  most  of  the  models  deal  with  specific  idealized 
geometries  and  packing  arrangements .  These  theories  provide  relations  between 
any  two  independent  elastic  constants  among  a  shear  modulus  (G),  a  bulk  modulus 
(K)  and  a  Young's  modulus  (E ) : 


K  =  K(Gm,  Km,  G,.  Kif  c) 
G  =  G(Gm’  V  V  V  c> 


wnere  subindices  m  and  i  refer  to  the  matrix  and  the  inclusion,  and  c  is  the 
volume  fraction  of  the  inclusion.  Smith  [11]  noted  that  through  algebraic 


manipulation,  most  of  these  tneones  can  oe  displayed  in  the  form  of  a  genera¬ 
lized  "rule  of  mixtures": 


where  G  and  j  are  the  shear  modulus  and  Poisson's  ratio  of  the  composite. 

All  of  the  previous  models  developed  by  elastic  analysis,  i.e.,  deformation 
of  the  matrix  and  inclusion,  have  a  common  feature:  a  specific  idealized 
geometry  and  packing  arrangement  has  to  be  assumed.  In  many  composite 
materials,  this  regularity  in  geometry  does  not  exist,  and  the  materia]  cannot 
be  identified  with  any  particular  arrangement  or  ensemble  valid  throughout  the 
whole  structure.  Thus,  the  previous  treatments,  in  the  strictest  sense,  do  not 
truly  repre;._nt  a  real  composite.  One  way  to  avoid  this  idealization  problem  is 


through  the  use  of  variational  techniques,  which  do  not  need  to  use  a  precise 
geometric  model. 

The  most  elementary  type  of  bounds  can  be  calculated  through  the  formalism 
of  the  complementary  energy  and  strain  energy  theorems  (Voigt  [12],  Paul  [13]). 
However,  the  resulting  bounds  are  often  too  far  apart  to  be  useful. 

In  order  to  improve  these  bounds,  Hashin  and  Shtrikman  [14]  formulated  new 
variational  theorems  on  the  basis  of  polarization  tensors.  However,  limita¬ 
tions  in  the  applicability  of  bounds  still  persist.  In  general,  it  is  not  known 
if  the  calculated  bounds  are  indeed  the  most  restrictive,  and  the  problem 
becomes  a  trade-off  between  tne  amount  of  information  introduced  in  the 
formulation,  and  the  resulting  contiguity  of  the  bounds. 

As  a  rule,  it  may  be  said  that  the  usefulness  of  bounds  is  restricted  to 
systems  where  the  properties  of  the  phases  are  similar,  i.e.,  metal  alloys.  In 
polymer  reinforcement,  however,  the  differences  between  matrix  and  inclusion  may 
be  consi derable.  As  a  consequence,  the  upper  and  lower  bounds  tend  to  separate, 
making  this  technique  inadequate  for  any  predictive  purpose. 

PROCEDURE 

As  indicated  earner,  a  basic  problem  in  the  theory  of  viscosity  of  suspen¬ 
sions  is  the  determination  of  correlations  useful  at  high  concentrations.  Most 
of  the  contributions  are  based  on  extensions  of  Einstein's  law  of  viscosity  at 
i nf i n1 te  di 1 uti on. 

Tne  analogous  problem  of  correlating  the  elastic  properties  of  heterogeneous 
systems  over  wide  ringes  of  concentration  has  not  been  successfully  solved.  As 
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will  be  discussed  later,  predictions  according  to  different  theories  show  a  con¬ 
siderable  lack  of  agreement. 

The  proposed  model  is  based  on  a  differential  scheme  similar  to  that  used 
by  Snnkman  [3]  and  Roscoe  [2]  for  the  prediction  of  the  viscosity  of  suspen¬ 
sions  at  high  concentrations.  In  a  later  publication,  Roscoe  [9]  tried  to 
extend  the  method  to  the  prediction  of  elastic  properties  of  composites. 
McLaugnlio  [30]  presented  a  similar  differential  scheme  and  evaluated  its  esti¬ 
mates  m  relation  to  the  Hashin-Shtrikman  bounds. 

Tne  present  model  is  based  on  tne  calculation  by  elastic  analysis  of  any 
two  independent  parameters  for  the  composite,  namely  the  shear  modulus,  the  bulk 
modulus,  ir  the  foung's  modulus.  Calculations  may  proceed  as  follows  [29]: 

1)  describe  the  deformation  state  (simple  shear,  etc.) 

21  choose  adequate  coordinates  for  a  given  geometry, 

3)  state  the  deformation  equations  for  matrix  and  inclusion, 

4)  state  continuity  of  stress  at  the  interface, 

5)  reduce  volume  integrals  to  surface  integrals  by  Eshelby’s 
formula  [5], 

6)  integrate. 

For  spnerical  inclusions  and  dilute  conditions,  the  shear  modulus  results  as 


15(1  -  yj  [1  -  (VGJ3 


Dvm  +  2(4  -5^)G1/G 
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c 


1  +  fc 


(3) 
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and  the  bulk  modulus  is: 


—  =  1  + 
K 


_ _ ( 1/^m)  (^i  ) _ 

1  *  C<Ki  -  +  WVV 


c  =  1  +  gc 


where  c  notes  the  volume  fraction,  and  the  subindices  m  and  i  refer  to  the 
matrix  and  the  inclusion. 

For  an  infinitely  rigid  inclusion  (G^  ®)  in  an  incompressible  matrix 

(v  =  0.5),  equation  (3)  results  in  a  relation  similar  to  Einstein's  law  of 
vi scosi ty : 


n  0 
m  m 


1  +  2.5  c 


Consider  a  suspension  containing  a  small  concentrati on  (volume)  c1  of  spheres 
into  which  is  placed  a  small  concentration  c 2  of  spheres,  and  into  this  again 
a  small  concentration  c^  of  spheres,  and  so  on,  up  to  n  sets  of  spheres,  and  let 
the  total  concentrati on  be  c.  Then  on  adding  another  set  of  spheres  of 
concentration  (volume  fraction)  cn+1>  the  new  total  concentration  (volume 
fraction)  is: 

c(1  '  cn+l}  +  cn+ 1  (6) 

So  the  increase  in  concentrat i on  is: 


Ac  cfinal  Cinitial  cn  +  l^ 


If  the  original  suspension  is  considered  a  homogeneous  medium  (matrix)  around 
the  new  spheres  which  are  in  small  concentrations,  it  is  possible  to  express 
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the  increment  in  shear  modulus  according  to  equation  (3)  as: 


Defining 


G  n+1 

n 


aG  =  Gn+1  '  Gn 


and  replacing  cn  +  ^  according  to  equation  (7), 


—  =  f  c  .  I  =  f— 
n  n  +  1  i 


'or  infinitely  small  increments  in  concentration. 


Ac  -*■  dc 


aG  *  dG 


dG  =  Gf  (9) 

d  c  ( 1  -c  ) 

This  differential  expression  for  the  increment  in  concentration  can  now  be 

used  in  equation  (3).  After  following  a  similar  reasoning  for  the  bulk  modulus 

K,  equations  13;  and  (4)  become: 

dG  =  G  1 5( i -v ) [  1-G ■ /G]  Q  . 

dc  [7-5v  +  2(4-5u  )G  ■  /G](  1  -c  ) 


~  K) 


^  . 1  *  [ ( ^ -j - k ) / ; !< * 1 4 / 3 ) g ) ] ;  (1-c) 


(11) 


where  v  - - —  is  the  roisson  s  ratio. 

2  ( 3  K  +  G ; 

Equations  (10)  and  (11)  constitute  a  coupled  system  to  be  numerically 
integrated  with  the  following  boundary  conditions: 


c  =  0 

G  =  Gm 

m 

(matrix  properties) 

K  =  K 

m 

(12) 

c  =  1 

G  =  G. 

i 

(inclusion  properties) 

K  »  K. 

1 

RESULTS 

1 '  Solid  Inclusions 

Smith  [11]  compared  the  behavior  of  well  known  models  for  the  prediction  of 
tie  mechanical  response  of  composites,  i.e..  Van  der  Poel  [lb],  Kerner  [17], 

3udi ansxy  [13],  etc.  In  order  to  test  the  performance  of  the  present  model,  the 
same  parameter  values  were  used  in  the  simulation.  Data  in  Figure  3a  have  been 
calculated  with  a  shear  modulus  ratio  G./G  =  30  and  Poisson's  ratios  v.  =  0.25 

i  m  i 

and  u  =  0.40,  representing  the  case  of  glass  spheres  embedded  in  an  epoxy 
ma  t  ri  <. 

Figure  3b  has  been  calculated  with  G./G  =  70,000,  v.  -  0.25,  ..  =  0.5, 

3  l  m  i  m 

representing  the  expected  shear  properties  of  glass  spheres  imbedded  in  a 
Iigntly  vulcanized  matrix  of  natural  rubber. 


In  both  cases,  a  considerable  discrepancy  among  the  different  models  is 


observed.  The  i nconvem ency  of  bounds  calculated  according  to  Hashin  and 
ShtrUman  L i a j  is  also  stressed.  The  predictions  of  the  present  simulation  are  close 
to  these  of  Van  der  Poel  [16]. 

Richard  (7)  presented  experimental  data  for  elastic  modulus  and  Poisson's 
ratio  of  a  polyester  matrix  reinforced  with  glass  spheres  (G^  =  5.94  10s  Pa, 

=  6.45,  3-  =  2.90  10iJ  Pa,  u-  =  0.21).  This  author  compared  his  experiments 
*;tn  the  predictions  of  several  theories.  Figure  4  indicates  the  excellent  pre¬ 
diction  of  tne  present  model,  while  tneories  by  Kerner  [17]  and  Hill  [20,21] 
only  perform  well  at  filler  concentrations  below  20  percent  volume. 

Richard's  data  [7]  were  also  used  by  Smith  [11]  for  comparison  with  other 

tneoretica1  contributions  as  indicated  in  Figure  5a  for  relative  moduli,  and 

Figure  5o  £or  Poisson's  ratio.  Calculations  were  performed  with  G^  =  5.84 

1 0 d  Ja,  v  =  0,44  G-  =  2.90  10‘°  and  v-  =  0.21. 
m  i  i 

Smitn  [20]  compared  his  own  experimental  findings  on  epoxy  resin-glass 


.meres  systems  with  several  theoretical  predictions.  The  properties  of  matrix 


; n  j  i  n c 1 u s 1  ir  are 


I  =  2.o3  1J9  Pa,  v  =  0.394,  E.  =  7.6  1010  Pa,  v.  =  0.23. 
m  m  i  i 


■  Its  those  predictions  are  plotted  and  compared  with  the  present  model 


o. 


-  go  at  ■  r1 


a^d  'll)  have  been  integrated  for  the  case  of  a  foam  by 


■  mole  mode  ,  for  me  incius’on. 
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Tne  simulation  was  compared  with  the  experimental  results  of  Gent  and  Thomas 
[24]  for  a  vulcanized  rubber  foam  (E^  =  2. 69  106  Pa,  =  0.49).  Figure  7  shows 
excellent  agreement  in  the  prediction  of  relative  modulus  versus  relative  den¬ 
sity  5, 

5  =  J-loam  =  1  -  c  (13) 

5  rubber 

where  c  is,  as  before,  the  volume  fraction  of  the  inclusion. 

The  performance  of  the  model  in  the  prediction  of  foam  properties  was  also 
tested  with  the  experimental  work  of  Moore  et  al.  [25]  for  a  variety  of  foamed 
thermoplastics  (PVC,  polypropylene  copolymer,  styrene-acrylonitri le  copolymer, 
etc.).  These  authors  observed  that  plots  of  relative  modulus  (relative  to  the 
solid  matrix)  versus  relative  densities  (6)  could  be  correlated  with  a  single 
curve  in  shear,  tension  and  compression  according  to: 


Figure  3  compares  the  predictions  of  several  theories  with  the  empirical 
findings  of  Moore  et  al.  [25]  and  the  results  of  the  present  model.  Simulations 
were  performed  for  a  polypropylene  copolymer  characteri zed  by  E^  =  1.13  GP'a, 

=  0.41.  Again,  the  proposed  method  shows  total  agreement  with  the  experimen¬ 
tal  values. 

3)  Low  modulus  incompressible  inclusions  (rubbers) 

For  the  limiting  case  of  an  incompressible  matrix  (v  =  0.5)  with  G-  =  0. 

m  i 


'foam  _  density  of  foam 

•solid  Ldensity  of  solid 


=  (6)2  =  (1  -  C ) 2 


equation  V1J)  admits  an  analytical  solution  as 


► 


which  represents  the  behavior  of  rubber  modified  polymers. 

Holliday  et  al.  [31]  compared  several  relevant  theories  for  a  case  repre¬ 
senting  soft  spheres  in  a  rigid  matrix  (Figure  9).  The  experimental  data 
correspond  to  acrylonitrile-butadiene-styrene  (ABS)  polymers  which  consist  of  a 
dispersion  of  soft  po lybutadiene  rubber  spheres  in  a  styrene-acrylonitri le  (SAN) 
copolymer  matrix  which  has  a  relatively  high  modulus. 

The  different  theories  examined  are  not  in  good  agreement.  The  proposed 
model  was  tested  with  typical  values  for  styrene-acrylonitri le  copolymer  (SAN) 
(Em  =  3.1  10'J  MPa,  =  0.38)  [32],  and  /Em  =  0.  Results  presented  in 
Figure  9  indicate  good  prediction  of  the  experiments. 

CONCLUSIONS 


Despite  many  contributions,  the  problem  of  predicting  the  elastic  behavior 
of  composite  materials  over  wide  ranges  of  concentration  has  not  been  satisfac¬ 
torily  solved. 

Different  theories  based  on  well  defined  geometric  models  present  a  con- 
sideraole  divergence  of  results.  The  use  of  bound  methods  has  been  proven  ina¬ 
dequate  fzr  problems  in  polymer  reinforcement,  where  solutions  become  a  trade¬ 
off  between  the  contiguity  of  the  bounds  and  the  amount  of  information 
'complexity)  introduced  in  the  model. 
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In  the  present  work,  a  correct  linearized  relation  for  the  elastic  defor¬ 
mation  of  matrix  and  inclusion  is  assumed  at  each  concentration.  This  relation, 
which  is  valid  locally,  can  then  be  extended  to  all  concentrations  under  the 
assumption  that  any  new  portion  of  filler  "sees"  the  existing  composite  as  a 
non-i nteracti ng  homogeneous  matrix. 

The  model,  presenting  no  adjustable  parameters,  accurately  represents 
available  experimental  data  on  solid  spherical  inclusions  and  foams.  The  method 
can  be  readily  extended  to  other  geometries  for  the  inclusion. 
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FIGURE  CAPTIONS 


Figure  1  Models  for  viscosity  of  suspensions. 

Figure  2  Models  for  prediction  of  elastic  response  in  composite  materials. 

Figure  3a  Relative  shear  modulus  versus  volume  fraction  of  filler  for  G^/G^  = 

30,  v,.  =  0.25,  *  0.4  (from  Smith  [11]).  Solid  curves,  model 

performance;  dashed  curves,  predictions  with  different  models: 

(1),  approximate  Van  der  Poel  [16];  (2),  Kerner  or  Hashin  and 

Shtrikman  highest  lower  bound  [14];  (3),  corrected  Van  der  Poel 

[16];  (4),  Budiansky  [18]. 

Figure  3b  Relative  shear  modulus  versus  volume  fraction  of  filler  for  G./G  = 

i  m 

70,000,  v.  =  0.25,  m  =  0.5  from  Smith  [11].  Solid  curve,  model 
performance;  dashed  curves,  predictions  with  different  models. 

(1),  approximate  Van  der  Poel  [16];  (2),  Hashin  and  Shtrikman 

highest  lower  bound  [14];  (3),  corrected  Van  der  Poel  [16]; 

(4),  Budiansky  [18]. 

Figure  4  Elastic  modulus  versus  volume  fraction  of  filler  for  the  case  of 

glass  beads  in  a  polyester  matrix.  Experimental  data  from  Richard  [19]. 
Solid  curve,  model  performance;  dashed  curves,  predictions  according 
to  1),  isostrain;  2),  isostress;  3),  Kerner  [17];  4),  Hill  [20,21]. 

Figure  5  Relative  Young's  modulus  E/Em  (5a)  and  Poisson's  ratio  (5b)  versus 

volume  fraction  of  filler  for  the  case  of  glass  beads  in  a  polyester 
matrix.  Experimental  data  from  Richard  [19].  Solid  lines,  model 
performance;  dashed  lines,  predictions  with  different  models. 

1),  Hashin  and  Shtrikman  highest  lower  bound  [14];  2),  Hashin  and 

Shtrikman  least  upper  bound  [14];  3),  Van  der  Poel  [16];  4), 


Budiansky  [18]. 


Figure  6  Relative  Young's  modulus  E/Em  versus  volume  fraction  of  filler  for 
the  case  of  glass  spheres  in  an  epoxi  matrix.  Experimental  data  by 
Smith  [23],  Solid  curve,  model  performance;  dashed  curves,  predic¬ 
tions  with  models  by  1),  Hashin  and  shtrikman  highest  lower  bound  [14] 
2),  Hashin  and  Shtrikman  least  upper  bound  [14];  3),  Van  der  Poel 

[16];  4),  Budiansky  [18]. 

Figure  7  Young's  modulus  of  a  highly  vulcanized  rubber  foam  E^  relative  to  the 
solid  rubber  E^  versus  volume  fraction  of  rubber  in  the  foam.  Solid 
line,  model  performance.  Experimental  points  were  determined  by  Gent 
et  al.  [24]. 

Figure  8  Relative  moduli  of  foams  versus  relative  densities.  Calculations 

here  have  been  performed  for  the  thermoplastics  described  by  Moore  et 
al.  [25].  Solid  curves,  model  performance  for  1),  Young's  modulus; 

2),  shear  modulus.  Dashed  curves,  predictions  according  to  different 
models.  3),  Kerner  [17];  4),  square-i n-square  model  [27];  5),  cube- 

in-cube  model  [27];  6),  law  of  mixtures.  Moore  et  al.  [25]  correlated 

experimental  moduli  with  a  square  law  which  results  coincident  with 
curve  1  (model  prediction). 

Figure  9  Relative  modulus  of  acrylonitri le-butadiene-styrene  (ABS)  copolymer 


versus  volume  fraction  of  polybutadiene.  Experimental  data  by 
Holliday  et  al.  [31],  Solid  curve,  model  performance.  Dashed  curves, 
predictions  with  different  models.  1),  Van  der  Poel  [16];  2)  Parallel; 
3),  Kerner  [17];  4),  Rei ner-Hashi n  [33];  5),  Series. 
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